Reentrant transitions in colloidal or dusty plasma bilayers 
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The phase diagram of crystaUine bilayers of particles interacting via a Yukawa potential is calcu- 
lated for arbitrary screening lengths and particle densities. Staggered rectangular, square, rhombic 
and triangular structures are found to be stable including a first-order transition between two dif- 
ferent rhombic structures. For varied screening length at fixed density, one of these rhombic phases 
exhibits both a single and even a double reentrant transition. Our predictions can be verified 
experimentally in strongly confined charged colloidal suspensions or dusty plasma bilayers. 
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Confined systems exhibit structural and dynamical be- 
havior very much different from the corresponding bulk 
state 0, y|- In particular freezing is strongly affected 
by the presence of a planar wall. In equilibrium, solid- 
ification near walls can occur at thermodynamic condi- 
tions where the bulk is still fluid (so called "prefreezing" ) 

Q . In nonequilibrium, the wall may act as a center 
of heterogeneous nucleation 5] in order to initiate crys- 
tal growth _6j. A system confined between two parallel 
planar walls exhibits various layered crystalline states at 
low temperature if the plate distance gets comparable to 
the mean interparticle distance. For hard spheres be- 
tween hard plates, geometric packing considerations lead 
to the stability of different crystalline lattices including 
multiple square and hexagonal layers 0] as well as buck- 
led [gl , rhombic H, |^ and prism superlattices . On 
the other hand, for pure coulombic systems like (classi- 
cal) electrons in quantum wells 11] or trapped ions 1121 . 
several crystalline bilayer-structures were reported |l3 l| . 

Most of our experimental knowledge of freezing in con- 
fining slit-like geometry is based on real-space measure- 
ments of mesoscopic model systems such as charged col- 
loidal suspensions between glass plates 0, 13 or of mul- 
tilayers of highly charged dusty plasmas [lj|. The ac- 
tual interaction between these mesoscopic "macroions" 
is neither hard-sphere like nor pure coulombic but is de- 
scribed by an intermediate screened Coulomb or Yukawa 
pair potential [isL due to the screening via addi- 
tional microions in the system. The screening length can 
be tailored by changing the microion concentration: for 
charged colloids, salt ions are conveniently added to the 
aqueous suspensions; the complex plasma, on the other 
hand, consists of electrons and impurity ions. 

In the present letter, we study the stability of differ- 
ent crystal lattices in bilayers of Yukawa particles as mo- 
tivated by the experimental model systems. The zero- 
temperature phase diagram is calculated for arbitrary 
screening lengths and particle densities We find a 

variety of different staggered solid lattices to be stable 
which are separated by either first or second-order phase 
transitions. The two known extreme limits of zero or 



infinite screening length corresponding to hard-spheres 
and the plasma 13, are recovered. For interme- 
diate screening lengths, the phase behavior is strikingly 
different from a simple interpolation between these two 
limits. First, there is a first-order coexistence between 
two different staggered rhombic lattices differing in their 
relative shift of the two unit cells. Second, one of these 
staggered rhombic phases exhibits a novel reentrant effect 
for fixed density and varied screening length. Depending 
on the density, the reentrant transition can proceed via a 
staggered square or a staggered triangular solid including 
even a double reentrant transition of the rhombic phase. 
All our theoretical predictions can in principle be verified 
in real-space studies of confined charged suspensions or 
dusty plasmas. 

In detail, our system consists of two layers containing 
in total N particles in the {x, y) plane. The total area 
density of the two layers is p = N/A with A denoting 
the system area in the (x, y) plane. The distance D be- 
tween the layers in the z direction is prescribed by the 
external potential confining the system. The particles 
are interacting via the Yukawa pair potential 



V{r) = 



exp(— Kr) 



(1) 



where r is the central separation. The inverse screen- 
ing length K, which governs the range of the interac- 
tion is given in terms of the micro-ion concentration by 
Debye-Hiickel screening theory. The energy amplitude 
Vo = Z'^ exp(2Ki?)K/e(l -I- scales with the square of 
the charges Z of the particles of physical hard core radius 
R reduced by the dielectric e permittivity of the sol- 
vent (e = 1 for the dusty plasma). Typically, Z is of the 
order of 100 — 100000 elementary charges such that V{r) 
at typical interparticle distances can be much larger than 
the thermal energy fc^T at room temperature justifying 
formally zero-temperature calculations. Then the energy 
scale is set by Vq alone and phase transitions in large bi- 
layer systems are completely determined by two dimen- 
sionless parameters, namely the reduced layer density 
?7 — pD^ /2 and the relative screening strength A = kD. 
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For zero temperature, the stable state is solid but differ- 
ent crystalline structures of the bilayers are conceivable. 
As possible candidate structures we assume that the two 
two-dimensional periodic lattices in the bilayers are the 
same, have a simple unit cell and are shifted relative to 
each other in the lateral direction by a displacement vec- 
tor c. If the two layers are labeled with A and B, the 
particle positions in the {x, y) plane of the two layers are 
given by 

Ryi(™,?^) = mai+na2, 

Rb(to,7i) = mai + na2 + c, (2) 

where ai and a2 are the primitive vectors of the two- 
dimensional lattice and m, n are integers. The total in- 
ternal energy U is obtained by the double lattice sum 

U^\Y1 V{\Ra-^'a\) + \ E V{\Rb-JMb\) 

+ E ViWB.A-'R-B^ + D^Y'^). (3) 

In the limit — s- oo, the stable crystalline structure min- 
imizes the total internal energy per particle U /N . 

We have minimized U /N with respect to ai, a2 and c 
under the constraint of prescribed density 77 for given A 
mapping out the phase diagram in the {rjX) plane. As a 
result, five typical staggered lattice structures turn out 
to minimize U /N for different 77. Adopting the notation 
developed for plasma bilayers [l8j, we label them by I, 
II, III, IV, and V. As summarized in Table HI phase I 
is the staggered rectangular crystal with a fixed aspect 
ratio 02/01 of -v/S; phase II is also staggered rectangular 
but with a different aspect ratio 7 interpolating contin- 
uously between phase I and the staggered-square phase 
III where 02/01 = 1. The staggered rhombic phase IV 
has two non-orthogonal lattice unit vectors (ai and a2) 
forming an angle 9 and contains a general lateral shift 
c = a(a.i + a2) between the two rhombic lattices. In 
fact, we find two possibilities for a defining two variants 
of stable rhombic phases which we call IVA and IVB. For 
IVA, a = 1/2 while a < 1/2 for IVB. Finally, phase V 



TABLE I: Structure and parameters of tlie different staggered 
bilayer crystals, ai is set to (ai,0) where ai is tlie nearest 
intralayer distance between particles. For phase II, 7 = a2/a\ 
is the aspect ratio. For phase IV, 9 is the angle between ai 
and sl2, and a is a free parameter characterizing the relative 
lateral interlattice shift c. 



Phase 


a2/ai 


c 


pal/2 


I. Rectangular 


(0, V3) 


(ai +a2)/2 


l/%/3 


II. Rectangular 


(0,7) 


(ai +a2)/2 


7 


III. Square 


(0,1) 


(ai +a2)/2 


1 


IV. Rhombic 


(cos 6, sin 6) 


(ai + 82)0: 


1/ sin 61 


V. Triangular 


(1/2,^3/2) 


(ai +a2)/3 


2/\/3 
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FIG. 1: Phase diagram of the Yukawa bilayer in the (r?, A) 
plane, (a) The hard sphere limit A ^ oo is sketched on top. 
The dashed (solid) lines denote continuous (discontinuous) 
transitions. The filled region corresponds to the coexistence 
domain of phases IV and V. The vertical arrow indicates the 
double reentrant behavior of phase IVA. The insets show the 
lattice geometries, where the filled (open) circles correspond 
to the lower (upper) layer, (b) Magnification of (a) show- 
ing a reentrant behavior of phase IVA occurring at moderate 
A. The four diamonds along the arrow indicate state points 
which were investigated by computer simulation at finite tem- 
peratures. 



is a staggered triangular crystal. Both phases III and V 
can be considered as special cases of the rhombic phase 
IV; the former has 6 = tt/2 and a = 1/2 while the latter 
is characterized hy — Tr/i and a — 1/3. 

The result for the phase diagram for a wide range of 
screening strengths (0 < A < 100; A — > oo) and densities 
(0 < rj < 0.8) is shown in Fig. ^ At v ery low screening 
A, we recover the known plasma limit [l8j, with our la- 
beling of the phases being in line with their sequence for 
increasing density rj. Phase I has a finite but extremely 
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small density region of stability up to 77 = 3.6 x 10~^ 
at A = [13 ■ For finite A, the I ^ II transition stays 
second-order and occurs at even smaller densities which 
decrease monotonically to zero as a function of A until 
the hard-sphere limit A — > cx) is reached. In this latter 
case ai is playing the role of an effective particle diam- 
eter. This is sketched by the vertical line in Fig. ^ 
The II — > III transition is second-order and the transi- 
tion densities increase drastically with growing A and in- 
terpolating monotonically between the plasma and hard 
sphere limit. More details of the I ^ II ^ III transition 
scenario are depicted in Fig. |2] where the aspect ratio 7 
of phase II is shown versus r] for different A. Phases I 
and III correspond to 7 = -s/S and 7 = 1, respectively. 
As can be clearly deduced from Fig. El the aspect ratio 
7 interpolates continuously as a function of rj between 
V3 and 1 for any A such that both the I ^ II and the 
II III transitions are second-order. In the hard sphere 
limit, 7 approaches 7/15 = —2rj + y^jf + S continuously 
which is also shown in Fig. [3 

Novel effects are observed for the III ^ IV — > V tran- 
sitions. First, for small A the III V transition proceeds 
via a IVA phase, the former being second and the latter 
first-order. For A « 8, however, there is a strong first- 
order transition directly from III to V with a large den- 
sity jump as determined by Maxwell's construction [20| . 
For even higher screening, A > 30, the III V transi- 
tion happens via the cascade III — > VIA — > IVB V. 
The stability range of the IVA phase becomes smaller for 
increasing A shrinking to zero in the hard sphere limit. 
Details of the III — > V transition scenario can be detected 
via the order parameters sm0 and a of the lattice min- 
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FIG. 2: Aspect ratio 7 — 02/01 for phase II versus density 
r] for different screening strengths A. The hard sphere case 
^hs is also shown. The lattice geometry is shown as an inset, 
where the filled (open) circles correspond to the lower (upper) 
layer. 
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FIG. 3: (a) Sine of the angle 9 and (b) the relative shift pa- 
rameter a = c/£ (with £ = |ai -l-a2|) versus density regarding 
the III ^ IV ^ V transition scenario for different A. The 
insets show the lattice geometry of phase IV. 



imizing the total potential energy at prescribed density 
rj. Plotting sm6 and a versus rj reveals the order of the 
transitions (sec Fig. O: A cusp, which is found for the 
III — + IVA transformation, implies a second-order tran- 
sition. All other transitions are first-order as signaled by 
discontinuous jumps in at least one of these order param- 
eters. The corresponding coexistence density gap is not 
shown in Fig. |21 but included in Fig. da). Across the 
IVA IVB transition the order parameter jump is small 
yielding a tiny density gap which can not be resolved in 
Fig. Uta). 

Our most striking result is reentrant behavior of the 
IVA phase at fixed density upon varying A as indicated 
in Fig. n]by the vertical arrow. For 0.5 < 77 < 0.525, 
there is reentrance of the VIA phase via the III phase. 
The full sequence over the whole range of A is IVA — > 
III ^ IVA IVB. For 0.530 < < 0.536 there is 
even a double reentrant behavior of the VIA phase via 



4 



the sequence IVA ^ V -> IVA -> III ^ IVA ^ IVB. 
This rich scenario is due to a subtle interplay of the 
range of the interaction in conjunction with the differ- 
ent bilayer lattice structures. Finally, at finite tempera- 
tures T, we performed extensive Monte-Carlo computer 
simulations with 800 particles in a rectangular-shaped 
box periodically-repeated in x and y direction and with 
hard walls of distance D in z-direction allowing fluctu- 
ating z-positions of the particles ( "buckling" ) . For fixed 
77 = 0.533, we investigated four states at A = 0.5, 3.0, 8.2, 
10 [see the four diamonds along the arrow in Fig. db)] for 
different T up to melting. The melting point is detected 
via a modified Lindemann-type criterion involving differ- 
ences of mean-square displacements of nearest-neighbors 
[1^ . We confirm that the reentrant behavior is stable 
with respect to increasing T up to melting. 

In summary, we have calculated the full phase dia- 
gram for a Yukawa bilayer at zero temperature by lat- 
tice sum minimizations. A competition between three 
length scales, namely the bilayer distance D, the aver- 
aged particle distance and the range 1/k of the 
interaction, induces a rich phase behavior which is dif- 
ferent from a simple interpolation of the extreme limits 
of the confined plasma and the hard sphere system. We 
predict a coexistence of two different rhombic phases at 
finite screening and a single and double reentrant scenario 
for one of the rhombic phases for varied "softness" of the 
interaction. These effects are in principle detectable in 
real-space experiments of charged colloidal suspensions 
confined between plates and in layers of dusty plasmas 
by tuning the screening strength via the microion concen- 
tration. The reentrant effect as obtained here in equilib- 
rium should also manifest itself as an interesting finger- 
print in nonequilibrium situations. For example, bilayer 
crystal nucleation and growth could be greatly stimu- 
lated via structures which are energetically close to the 
stable ones [l^l- Soft particle interactions different from 
the Yukawa type of Eq. (^), as e.g. inverse power poten- 
tials where V{r) oc r~", will lead to similar reentrant 
effects as long as the softness of the potential (e.g. the 
exponent n) is varied. Different realizations of soft in- 
teractions occur in sterically-stabilized colloids, in spher- 
ical block-copolymer micelles and in star polymers and 
dendrimers, where the softness of the effective interac- 
tion can be tuned by the length and grafting density of 
the polymer chains or the solvent quality 22]. Hence 
the reentrant scenario should also occur in foam films 
containing polymer bilayers |23j . Finally, for a general 
external potential confining the particles to layers, the 
bilayer distance D is not prescribed but the system will 
minimize its total energy realizing an optimal D. In this 
case, second-order phase transitions will still be described 
in terms of scaled parameters. This implies a universal 
behavior of our bilayer phase diagram. In a general ex- 



ternal potential, however, the system has the additional 
possibility to split into tri- and higher order multilayers. 
This can happen either discontinuously or continuous via 
merging prism phases. Details have to be explored in fu- 
ture studies. 
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